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1. Introduction 

The investigation of infrared properties of Landau-gauge gluon and ghost propagators — in 
order to test predictions of the Gribov-Zwanziger confinement scenario, or of the Kugo-Ojima 
scenario — has generated a flurry of papers in the last couple of years. While the formulation 
of these two scenarios, their equivalence (or lack thereof) and the chosen procedure for solving 
these propagators' Dyson-Schwinger equations have become matters of heated debate, a consis- 
tent picture emerges from lattice studies. Based on this picture, which is in partial disagreement 
with the predictions of both scenarios, present activity focuses on critical revisions of the origi- 
nal scenarios, on discussion of their main assumptions and/or implications, and on whether or not 
there are physical criteria to prefer a solution of the Dyson-Schwinger equations (or other func- 
tional methods) of the "scaling" or of the "massive" type. In this way, one hopes to have gained a 
deeper insight into confinement in Landau gauge. We will not review these various analytic studies 
here, but rather refer to recent status reports and overviews contained in these proceedings, such as 
Ref. [0], which summarizes the so-called refined Gribov-Zwanziger framework, Refs. [§] and [g], 
which address the problem of the characterization of the Kugo-Ojima scenario and its relation to 
the Gribov-Zwanziger one, or Ref. which reviews the solution of Dyson-Schwinger equations 
in Landau gauge and discusses phenomenological applications. We do, however, attempt to review 
thoroughly the recent literature on lattice studies of the topic. 

In what follows we consider tests of the original Gribov-Zwanziger confinement scenario, i.e. a 
vanishing gluon propagator and an enhanced ghost propagator in the infrared limit. On the contrary, 
in the refined Gribov-Zwanziger framework mentioned above, one expects a finite (nonzero) gluon 
propagator and a free ghost propagator in the same limit. These two cases are consistent respec- 
tively with the so-called scaling and massive (or decoupling) solutions of the Dyson-Schwinger 
equations. The latter behavior has been strongly favored in all recent extensive lattice studies, 
performed on very large lattices for pure SU(2) and SU(3) gauge theory. As argued in []|], the es- 
sential features of the original Gribov-Zwanziger confinement scenario are not incompatible with 
these findings, since violation of reflection positivity is clearly observed for the gluon propagator 
and enhancement with respect to the free propagator is seen for the ghost at intermediate momenta. 
Also, a logarithmic enhancement of the ghost propagator might be possible in the continuum limit. 
A scaling solution (with nontrivial infrared exponents), however, is ruled out. 

From the point of view of lattice simulations, we must strive to keep under control the various 
sources of systematic errors that might obscure the true infrared behavior of the propagators, in or- 
der to conclude that the behavior described above is firmly established. Of course, once the massive 
behavior is confirmed from the simulations, one must understand why this behavior arises in the 
infrared limit of the theory, and how we might reconcile it with a confinement mechanism. Again, 
insight into the problem can hopefully come from the numerical simulations themselves. Also, 
it should be noted that the comparisons mentioned above assume that the gauge definition on the 
lattice is physically equivalent to the continuum one, an issue that should be carefully investigated. 
We will refer here to the so-called minimal Landau gauge condition 

In Section |2] we review some aspects of the numerical simulations, with special attention to the 
main possible sources of systematic errors. We also attempt an overview of references and chronol- 
ogy of the main recent lattice results on the topic. In Section || we summarize interesting constraints 
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on the infrared behavior of the propagators, written in the form of upper and lower bounds at fixed 
lattice volume. Although these bounds were introduced as a guide to the infinite-volume extrap- 
olation, we believe they will be useful tools to investigate why the propagators have the observed 
behavior, since they naturally relate to a statistical interpretation of the gluon propagator and to 
a clearer view of the ghost propagator in terms of the spectrum of the Faddeev-Popov operator. 
Section |] is dedicated to the results from very large lattices mentioned above. We summarize the 
analysis of data from our simulations of the pure-SU(2) case, which are essentially equivalent to 
the corresponding results by other lattice groups. In Section || we address the /3 = case, where 
various sources of systematic errors may be investigated more easily. The observed behavior is 
compared with the one at finite j3. Finally, we present our conclusions in Section [6[ 

2. The simulations 

On the lattice, the gauge action is written in terms of oriented plaquettes, formed by the link 
variables U^{x), which are elements of the gauge group SU (N c ) and change under gauge transfor- 
mations as 

C/ M (x) -»• [/*(*) = g(x)U^{x)g(x+fl)\ (2.1) 

where g G SU(N C ), lattice sites are labeled by x and pL refers to the directions along the lattice. 
Consequently, all closed loops are gauge-invariant quantities, including the plaquettes in the pure- 
gauge-theory action 

S({U}) = £ £ £9iTr[U- U^(x)U v (x+fl)U;(x+v)ul(x)) . (2.2) 

Here /3 is the lattice parameter, related to the bare coupling constant as j8 = 2N c /g^. Except 
where otherwise indicated, we consider (symmetric) hypercubic lattices of four space-time dimen- 
sions. The lattice size in physical units is given by L = Na, where N is the number of points per 
lattice direction and the lattice spacing a is expressed in physical units. The physical volume is 
thus V = L 4 . (Note that one often refers to the "lattice volume" N 4 = V /a 4 .) We assume periodic 
boundary conditions. 

Let us remark that the resulting path integral has a finite (group) integration volume and there 
is in principle no need for gauge fixing on the lattice. A procedure for fixing the gauge numerically 
may nevertheless be implemented in a straightforward way in the simulation, without the need to 
consider the Faddeev-Popov matrix For minimal Landau gauge, one minimizes the functional 

= 1 " z^II 9mf W (2 - 3) 

with respect to the gauge transformations {g}. Indeed, the first derivative of this functional yields 
the familiar Landau gauge condition of null (lattice) divergence of the gauge field, which is defined 
in terms of the link variables as 

M*) = ^[^W-^W] traceless" ^ 
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The prescription for fixing the gauge in a simulation is thus to: 1) produce a gauge-link configura- 
tion {U } as usual, 2) to find {g} that is a (local) minimum of the functional $ (holding {U } fixed) 



and 3) to transform {U} following Eq. ( |2.1[ ) for the selected {g}. The resulting configuration is a 
gauge-fixed realization of the link variables, which will be used to compute observables of interest 
such as the gluon and ghost propagators. 

The gluon propagator is given in Landau gauge simply by 

Ofvip) = Z e ~ 2i * P ' X ( A » (*)^(0)> = S ab (g*v - Ej ^ L ) D(P 2 ) , (2-5) 

where p is the momentum and a, b are color indices. It is therefore determined solely by the scalar 
function D(q 2 ) associated to its transverse component. In the original Gribov-Zwanziger scenario, 
gluon confinement is associated with violation of reflection positivity for the gluon propagator (in 
real space-time). 

As pointed out above, lattice gauge fixing is accomplished without the need to compute the 
Faddeev-Popov matrix „#. Nevertheless, the matrix can be obtained directly from the second 
variation of the gauge-fixing functional S, which corresponds to the Jacobian of the gauge-fixing 
condition. It is interesting to note that in this way there is also no need to consider the ghost fields 
explicitly [Q]. The ghost propagator G(p 2 ) is given by the inverse of M as 

i -lnip\x-y) 
G{ P ^ = W~i I v (a,x;a,y)) . (2.6) 

™c 1 x,y,a ' 

An infrared enhancement of G(p 2 ) with respect to the tree-level ghost propagator G(p 2 ) ~ p~ 2 
is expected in the original Gribov-Zwanziger scenario (and in the Kugo-Ojima one) as a sign of 
confinement. 

Let us note that the known problem of Gribov copies is present on the lattice as well, since 
each local minimum of the functional $ corresponds to an equivalent (lattice) gauge copy. The 
algorithm for fixing {g} has in principle no control over which copy gets selected. Because of the 
minimization, we know that ^# is positive semi-definite and, as a result, the sampled copies are 
inside the first Gribov horizon Q, which is delimited by the vanishing of A m ,„, the smallest nontrivial 
eigenvalue of ^ . It is usually argued that a unique copy might be obtained, corresponding to the 
fundamental modular region A, if one were able to determine the global minimum of S. This region 
has been studied on the lattice in [S]. In any case, Q. and A are shown to be convex regions of very 
high dimensionality, which likely constrains the statistical weight of gauge configurations to lie 
near their boundary. In particular, one should check if the sampled configurations have vanishing 
Xmin as the lattice volume goes to infinity. 

Lattice simulations have been carried out since the mid 1980s for the gluon propagator [^] 
and since the mid 1990s for the ghost [10] (see also [[TT|]). Early studies have established that the 



gluon propagator is not enhanced at small momenta Jl2|], but did not allow further conclusions 
about its infrared behavior. A turnover point in momentum (suggesting a null propagator) could 



only be seen at strong coupling [13] or in three space-time dimensions [14]. However, even in 
this simplified case, a later study on a very large lattice (of volume 140 3 ) was still not conclusive 



[15], although it was possible to fit D(0) to zero in an infinite- volume extrapolation, and violation 



of reflection positivity was clearly seen [|T6|]. For the ghost propagator, infrared enhancement was 
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observed, but the corresponding infrared exponent seemed to become smaller as lower momenta 
became available ^ [TlJ, |l7[ [if]. It was also clearly shown that A mi „ goes to zero with increasing 



lattice volume [ 19| , 18 ]. These studies were complemented by investigations of the strong coupling 



constant (see e.g. Q20| , j21fl ) and several three-point vertices (see e.g. [ |22| , [23[ [24] , |18|]). This was the 
status until 2006. As became clear later, systematic effects had not yet been properly taken into 
account, which limited the conclusions (or the lack thereof) of these studies. We now pause for a 
moment and list the main such possible effects (and related references) below. 

• Gribov-copy effects: this is a very important issue. As commented above, usual simulations 
do not take fluctuations in the values of the propagators due to Gribov-copy effects into 
account. A few studies have considered the determination of the absolute minimum of the 



gauge-fixing functional (see e.g. Qll| , [25|]) or other criteria to fix the gauge [g6J. It was 
generally found that the effect of Gribov copies decreases as the lattice volume increases. 
This statement must be taken with a grain of salt, since the number of copies surveyed is 
limited and one does not know for sure if this number is large enough at a given volume 
to allow the determination of the global minimum. Studies of the exact structure of Gribov 
copies are now being carried out (on small lattices) [|27|]. Let us mention that it was argued by 



Zwanziger [ |28| ] that averages taken in the fundamental modular region should coincide with 
averages in Q. in the infinite- volume limit. We thus conclude that a sign of significant Gribov- 
copy effects has not yet been seen and the effects observed so far are probably connected to 
the next item below. (We do note, however, that a very recent study has reported on sizeable 



effects at large lattice volumes [29].) 



Finite-volume effects: perhaps surprisingly, these are the most serious systematic effects we 
have to deal with. To be sure, lattice simulations must be carried out at finite lattice volumes, 
since computers have finite memory. As mentioned above, the physical extent of the lattice L 
is given by the number of lattice points along each direction multiplied by the lattice spacing 
a in physical units, which is directly related to the lattice parameter j3. To simulate closer 
to the continuum limit one must go to smaller a, or equivalently to larger j8, while keeping 
the simulated lattice large enough to represent the relevant energy scales of the problem. 
Strictly speaking, one would need an extrapolation to infinite lattice volume at each fixed 
value of /3 . (A continuum extrapolation would additionally require running at increasingly 
smaller values of a.) In usual lattice applications, though, taking the infinite-volume limit is 
not among the most serious issues, since one typically just needs to have a sufficient number 
of points to ensure a physical lattice size of the order of the relevant hadronic scale, i.e. 
around 1 fm. The main effort is then to go to very small a, in order to avoid discretization 
errors (addressed below). In studies of the infrared limit, however, the situation is different, 
and finite-size effects play an important role. This happens because the infrared limit lies at 
small p, corresponding to large L. (Note that the smallest nonzero momentum that can be 
represented on a lattice of side L is ~ 2n/L.) 

Discretization effects: as indicated above, the effects due to simulating at nonzero lattice 
spacing a are not expected to be so serious in the infrared limit, because the energy scale 
associated with the cutoff a (which is ~ I /a), is sufficiently high compared to the typical 
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momenta of interest. In other words, the long wavelengths we are interested in do not re- 
solve the lattice spacing and are not much affected by it. Nevertheless, discretization errors 
may be important for the breaking of rotational symmetry as well as for the possible different 
discretizations of the gluon field and of the gauge-fixing condition. There are ways to reduce 
effects due to the breaking of rotational symmetry, such as cutting out the momenta charac- 



terized by large effects [12] (the so-called cylindrical cut), improving the lattice definition of 



the momenta [30] and including (hypercubic) corrections into the momentum-dependence of 



the Green's functions [31]. As for the discretization of the gluon field and of the lattice Lan- 



dau gauge condition, several different definitions may be considered (see e.g. [32]). These 
studies have usually found that different discretization procedures lead to gluon propagators 
that differ only by a multiplicative constant, which can be reabsorbed in the (multiplicative) 
renormalization of the propagator. 

Unquenching: in the Gribov-Zwanziger and related scenarios, one hopes to get an under- 
standing of confinement in the static-quark limit, where there is no string breaking and the 
confinement problem may be phrased as a search for explaining why an area law develops 



[]33[]. Thus, it should be sufficient to consider the pure-gauge theory treated here, also known 
as the quenched approximation. Nevertheless, an important question is how the picture gets 
affected once dynamical quarks are introduced in the simulations. Studies done so far (on 
relatively small lattices) show qualitatively the same behavior as in the pure-gauge case [B4p. 



In 2007, studies of Dyson-Schwinger equations on the torus [ ]35[ ] hinted that finite-size effects 
might indeed be plaguing results from lattice simulations and predicted that physical lattice sides 
of the order of 15 fm might be needed to begin to see the expected (conformal scaling) infrared 
behavior of the propagators. At about the same time, two other predictions of similar studies 
were verified in simulations: the (quantitative) equivalence of infrared propagators in the SU(2) 



and SU(3) cases 1 [37, 38] and the verification of conformal scaling behavior in two space-time 



dimensions [p9Q. That same year, three groups came out with studies on very large lattices, which 
were all presented at the Lattice 2007 conference. The Berlin-Dubna group considered 80 4 lattices, 
corresponding to a lattice extent of 13 fm, in the SU(3) case (Their study of the gluon 



propagator was later extended to 96 4 lattices, corresponding to a lattice extent of 16 fm [41].) The 
Adelaide group considered 112 4 lattices, corresponding to 19 fm, in the SU(2) case [pq]. Similarly, 
we considered 128 4 lattices (corresponding to 27 fm) in the SU(2) case, plus three-dimensional 



lattices of size 320 3 , corresponding to 85 fm []42[]. What these studies showed was puzzling. On 



the one hand, the large volumes clearly allowed a better view of the infrared picture; on the other, 
this view was nothing like what the authors had imagined it would be! In fact, going to large 
volumes in the hopes of seeing a null infrared gluon propagator not only established that D(0) was 
not null (even in three dimensions), but also exposed the fact that the previously seen enhancement 
of the ghost propagator goes away at very small momenta, and the data are consistent with a flat 
ghost dressing function in the infrared limit. 



'Although one cannot expect this to hold at high values of p, it is conjectured that the cases N c = 2 and N c = 3 
have the same infrared behavior. A recent comparison along a wider range of momenta presented in [^] shows some 
discrepancies between the two cases. 
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This behavior is analyzed in Section ^ In the next section, on the other hand, we comment on 
rigorous bounds for the propagators, introduced as a guide to the infinite- volume extrapolation. 

3. Bounds on propagators and statistical interpretation 

As discussed in the previous section, one of the main difficulties in the lattice simulations is 
the extrapolation of gluon- and ghost-propagator data to infinite lattice volume. In fact, the correct 
volume dependence of the data may not be easily inferred from the behavior on medium-size (or 
even very large) lattices, especially since some quantities, such as the zero-momentum gluon prop- 
agator, are quite noisy. It is then very helpful to obtain constraints on the infrared behavior of the 
propagators, as the upper and lower bounds discussed in this section. We remark that these bounds 
are valid at each lattice volume and must be extrapolated to infinite volume, just as for the propa- 
gators. The advantage is that the bounds are written in terms of quantities that are more intuitive 
than the propagators themselves, making it easier to guess the expected volume dependence of the 
propagators and possibly allowing an explanation of the infrared behavior observed in the data. 



In the case of the gluon propagator, we obtain [ [431 ] the bounds 



V (M(0)} 2 < D(0) < Vd{N 2 c - 1) (M(0) 2 ) , (3.1) 



where d is the dimension, V is the volume, 



D(0) = ;?(]vFT)M (0)|2) M 



is the zero-momentum propagator and M(0) is defined as 

1 

V ■ ' h.u 

Let us also define the "magnetization" 

1 



From the above definitions we can get a statistical interpretation for the quantity on the right-hand 



side of (3.1): it is essentially the susceptibility associated with the magnetization M'(0) (since the 



average of this magnetization vanishes, due to the residual global gauge symmetry). By analogy 
with a <i-dimensional spin system one would thus expect to see V (M(0) 2 } ~ const, i.e. the statistical 
variance of the magnetization is proportional to the inverse of the volume, a behavior known as self- 
averaging. At the same time, considering the statistical fluctuations in the Monte Carlo sampling 
of M(0), we would expect (M(0)} 2 to have the same volume dependence as (M(0) 2 } The 
simple statistical argument presented above suggests that both (M(0)) 2 and (M(0) 2 ) should show 
a volume dependence as 1/V. On the other hand, this suppression with l/V is compensated by 



the volume factor for both bounds in Eq. (3.1). Consequently, if this suggested behavior for the 
susceptibilities is verified, D(0) converges to a nonzero constant in the infinite- volume limit. As 
explained in the next section, this is what one observes in the simulations. Note that the bounds 
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in Eq. (3.1) apply to any gauge and that they can be immediately extended to the case D(p • ) with 



p ^ 0. We also note that a very interesting stochastic interpretation for the gluon propagator has 



been investigated in [440. 



Also in the case of the ghost propagator, a more intuitive picture comes from noticing that in 



Landau gauge, for any nonzero momentum p, one finds [ ]45[ ] 

iJ—t^-IslYmin^p)? < G(p 2 ) < (3.5) 

iV c 1 A lmn a A m j n 

where X min is the smallest nonzero eigenvalue of the Faddeev-Popov operator jtft and \jr min (a, p) is 
the corresponding eigenvector. Note that the upper bound is independent of the momentum p. If 
we now assume A mi „ ~ L~ v and G{p 2 ) ~ p~ 2 ~ 2K at small p, we have that 2 + 2k < v, i.e. v > 2, 
is a necessary condition for the infrared enhancement of G(p 2 ). A similar analysis can be carried 
out for a generic gauge condition. Consider the Gribov region Q, where all eigenvalues of ^ are 
positive. In the infinite-volume limit, as mentioned above, entropy favors configurations near the 



Gribov horizon dQ, where A m , n goes to zero. Thus, inequalities such as ( |3.5| ) can tell us if one 
should expect an enhanced G(p 2 ) when the Boltzmann weight gets concentrated on dQ. (This 
answers the question posed in p6|].) In other words, it is clearly necessary — but not sufficient — 
to have a vanishing Ami„ as the volume tends to infinity in order to observe enhancement of G{p 2 ). 
The upper bound in Eq. ( |3~3| ) was tested for our data in p5|], albeit with limited statistics for A m ;„. 
We find V « 2, consistent with finding K — > from fits of G(p 2 ), as shown in the next section. 

Let us also mention that a possible connection between the infrared behavior of the gluon 
propagator and the appearance of nontrivial zero modes of ^# for configurations near dQ has been 



recently presented in [47]. The two types of scenarios obtained there can probably be related to the 



massive and the scaling solution for the propagators. 
4. Huge lattices 

Our study in the SU(2) Landau case [|3|], using the very large lattices mentioned above, is 
summarized here. In Fig. [T] we show data for the gluon propagator in d = 4 and d = 3 for a wide 
range of (large) lattice volumes, indicating that D(0) remains nonzero in the infinite- volume limit. 
Similar results are obtained in [ ]3~8[ , |4C| , BT| ], as mentioned before, but also in [48], which takes 
Gribov-copy effects into account and in [09], which uses improved actions and anisotropic lattices. 
We have investigated the volume dependence of the bounds in Eq. (|3~l|) and found remarkably 



good agreement with the predicted 1/V behavior for (M(0)) 2 and (M(0) 2 ), thus implying a finite 
nonzero value for D(0) in the infinite-volume limit. More precisely, by fitting the two quantities 
to l/V a we get the exponents a respectively 0.995(10) and 0.998(10). A similar analysis for the 
SU(3) case (considering somewhat smaller volumes) yields the exponents 1.058(6) and 1.056(6) 



[50]. Violation of reflection positivity for D(x) is seen in all cases. 

Our data for the ghost propagator support a tree-level (or free) form in the infrared limit. This 
behavior is better seen if one considers the dressing function p 2 G(p 2 ), as shown in Fig. ||. Indeed, 
the data can be well fitted [^] by the form a — b [log ( 1 + cp 2 ) + dp 2 ] / ' ( 1 + p 2 ) , consistent with 
K = in the infrared limit. Note also that, for smaller p, this form is equivalent to a — Mog(l +cp 



2a 



(proposed in [51]), where c may be related to a gluon mass. This is also observed in d = 3. We 
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0.5 1 1.5 2 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



P P 

Figure 1: The gluon propagator D(p 2 ) as a function of the lattice momenta p (both in physical units) for 
the pure-5t/(2) case in d — 4 (left), considering volumes of up to 128 4 (lattice extent ~ 27 fm), and d = 3 
(right), considering volumes of up to 320 3 (lattice extent ~ 85 fm). 



remark again that enhancement is seen at intermediate momenta and that, depending on how the 
fit parameters change with the lattice spacing, there might be a logarithmic enhancement in the 
continuum limit. 



5. p = 

An interesting laboratory to test for various sources of systematic errors in the simulations is 
the apparently trivial case of /3 = 0, i.e. no dynamics from the lattice Boltzmann weight associated 



with the action [see Eq. (p.lQ]. This was considered in [13, 52, 53]. Since j3 = corresponds to 



an unphysical limit, the issue of setting the lattice scale (given by a) is delicate. As discussed in 



detail in [53], a possible choice is a — > °°. This is convenient because the lattice extent will already 
be infinite and there will be no finite-size effects. On the other hand, if a is large and p is not very 
small there could be discretization effects. We must see from the data which effect is predominant. 
(The latter effect may also be measured in terms of breaking of rotational invariance.) As shown 
in [^3|], we see: 1) clear violation of reflection positivity for the gluon propagator, 2) a seemingly 



finite and nonzero limit for D(0), including analysis with the bounds in Eq. (3J_), 3) a very good fit 
of p 2 G(p 2 ) to the form a — Mog(l + cp 2 ) and 4) no finite-volume or rotational-symmetry-breaking 
effects, suggesting that the data are in the deep infrared limit and at infinite volume. Essentially, 
we see the same infrared behavior as for finite j3. Note that our conclusions differ somewhat from 



[52] 



6. Conclusion 

Lattice simulations of infrared Landau gluon and ghost propagators have come a long way 
in the past couple of years. The current paradigm is that of a massive gluon and a free ghost, as 
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4.5 I — . ■ — h — . 1 h — . 1 n — i ■ h — . ■ n 5 




1 1 — 1 — 1 — 1 — 1 "J — 1 1 



0.001 0.01 0.1 1 10 100 0.0001 0.001 0.01 0.1 1 10 

P P 

Figure 2: The ghost dressing function p 2 G(p 2 ) as a function of the lattice momenta p 2 (both in physical 
units) for the pme-SU (2) case in d — 4 (left), considering volumes of up to 128 4 (lattice extent ~ 27 fm), 
and d = 3 (right), considering volumes of up to 320 3 (lattice extent ~ 85 fm). Note the logarithmic scale for 
p 2 . We also show the fits to a — b [log(l +cp 2 ) +dp 2 ~\ /(l + p 2 ). For d = 4, the fit parameters a,b,c,d are 
respectively 4.32(2), 0.38(1), 80(10), 8.2(3). 



proposed long ago by some Dyson-Schwinger-equation studies (see e.g. [54]). These simulations 
are greatly influenced by interchange with researchers who use analytic and semi-analytic methods, 
and vice-versa. This synergy has contributed to making infrared QCD a very productive field, with 
an active community and several dedicated workshops every year. 
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